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We calculated the effects of spin-orbit interaction (SOI) on the energy bands, ballistic conductance 
(G) and the electron-diffusion thermoelectric power (^d) of a nanowire by varying the temperature, 
electron density and width of the wire. The potential barriers at the edges of the wire are assumed 
to be very high. A consequence of the boundary conditions used in this model is determined by 
the energy band structure, resulting in wider plateaus when the electron density is increased due 
to larger energy-level separation as the higher subbands are occupied by electrons. The nonlinear 
dependence of the transverse confinement on position with respect to the well center excludes the 
"pole- like feature" in G which is obtained when a harmonic potential is employed for confinement. 
At low temperature, Sd increases linearly with T but deviates from the linear behavior for large 
values of T. 

PACS numbers: 73. 23. Ad, 75.70.Tj,81.07.Gf 

I. INTRODUCTION 

The electronic transport and photonic properties of a two-dimensional electron gas (2DEG) such as that found at 
a semiconductor heterojunction of GaAs/AlGaAs have been the subject of interest and discussion for many years 
now.i Related physical properties of narrow quantum wires of 2DEG have also been the subject of experimental 
and theoretical investigations because of their potential for device applications in the field of nanotechnology.— It 
is thus necessary to specify the model for the edge of a narrow quantum wire.— ^— Here, we analyze the role played 
by the boundaries on the ballistic electron transport in a nanowire of 2DEG where the Rashba spin-orbit interaction 
(SOI) is included. The role of SOI on collective properties of the 2DEG has been investigated.—"— The quasi-one- 
dimensional channel may be formed by applying a negative bias to a metal gate placed on the surface and depleting 
the 2DEG below it. Alternatively, the quasi-one dimensional channel may be made by etching all the way down to the 
active layer. In either way, it has been demonstrated that the width of the wire could influence the electron transport 
properties. In the latter case, the edges are sharper and may be appropriately modeled by sharp and high boundary 
conditions. More recently, Brey and Fertig— have pointed out that graphene nanoribbons have different collective 
plasmon dispersion relations depending on whether they have armchair or zigzag edges. 

It is well established that the spin-orbit coupling is an essentially relativistic effect: an electron moving in an 
external electric field sees a magnetic field in its rest frame. In a semiconductor, the interaction causes an electron's 
spin to precess as it moves through the material, which is the basis of various proposed "spintronic" devices. In 
nano-structures, quantum confinement can change the symmetry of the spin-orbit interaction. The relativistic motion 
of an electron is described by a Dirac equation. These effects combine to form both an electric dipole moment and 
the Thomas precession which is due to the rotational kinetic energy in the electric field.— The two mechanisms 
accidentally have very close mathematical form and consequently combine in a very elegant way. The SOI Hamiltonian 
can be obtained from the Dirac equation by taking the non-relativistic limit up to terms quadratic in v/c. This limit 
can be achieved either by expanding the Dirac equation in powers of v /c or by making use of the asymptotically exact 
Foldy-Wouthousen transformation. 

We include the effects due to edges through sharp and high potentials at the boundaries. As a result, we are not 
able to solve the Rashba SOI model Hamiltonian to obtain analytic solutions for the eigenenergies and eigenfunctions. 
The reason for this is due to the fact that the solution manifestly contains quantum interference effects from multiple 
scattering off the edges. We solved the eigenvalue problem numerically, obtaining the energies as a function of the wave 
vector ky parallel to the edge of the nanowire shown schematically in Fig.[T] This model for the edges is different from 
that employed in previous works.—"— We calculate the ballistic conductance and the electron-diffusion thermoelectric 
power for this quasi-one-dimensional structure by assuming that the length of the channel between source and drain 
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is less than the electron mean free path. In addition, we assume that the width of the channel is of the order of the 
de Broglie wavelength. 

The outline of the rest of this paper is as follows. In SecHIl we present our model for the wide quantum wire with 
spin-orbit coupling and specified boundary conditions to simulate the effects arising from the edges of the nanowire. 
We also present numerical results for the energy bands in order to study the combined effect of the boundaries and 
the SOI. Section IIIII is devoted to a brief description of the way in which our calculations are done for the ballistic 
conductance and electron-diffusion thermoelectric power when the energy bands are symmetric with respect to the 
wave vector ky parallel to the edges of the nanowire.— Numerical results and discussion of the conductance and the 
electron-diffusion thermoelectric power as functions of electron densityand temperature, for various wire widths and 
Rashba parameters, are given in Sec. lIVI A summary of our results is presented in SecfVl 

II. MODEL FOR ENERGY BAND STRUCTURE 

It is now well established that, spin-orbit coupling is an essentially relativistic effect. The relativistic motion of 
an electron is described by the Dirac equation that contains both effects (electric dipole and Thomas precession) in 
the spin-orbit interaction and does so in a very elegant way (see, e.g., textbooks— i^i). The SOI Hamiltonian can be 
obtained from the Dirac equation by taking the non-relativistic limit of the Dirac equation up to terms quadratic in 
v/c inclusive. This limit can be attained in two different ways: by direct expansion of the Dirac equation in powers 
of v/c and by the asymptotically exact Foldy-Wouthuysen transformation—. The Hamiltonian for an electron in the 
quadratic [0(w^/c^)] approximation is the sum 

n = nso + ^H, (1) 

where AH is the free-particle Hamiltonian and 
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describes the SOI within the material and includes both contributions to the spin-orbit coupling from the electric 
dipole and the Thomas precession (caused by the electric field) mechanisms. This result is general since it was derived 
from the Dirac equation, an exact relativistic equation for the electron, and includes all possible relativistic effects, 
whatever might be their kinetic source. When an electron gas at a heterojunction is confined to the xy-plane so that 
the electrostatic potentialis spatially-uniform along the heterostructure interface and varies only along the z axis, the 
Hamiltonian in Eq. contains just the contribution arising from its confinement along the z direction. For a quasi- 
one-dimensional structure, a second term must now be added to account for the extra local confinement produced by 
the electric field within the xj/-plane. 

For quantum wires, the width of the potential well is comparable with the spatial spread of the electron wave 
functions in the z direction. Therefore, in order to determine an effective electric field acting on electrons in the 
potential well, one should calculate an average of the electric field E{z) over the range of the z variable where the 
wave function is essentially finite. Consequently, one can model the averaged electric field by a potential profile. 
In principle, all potential profiles can be classified in two ways. In the first case, the average of E{z) is negligible 
although E{z) itself may not be zero or even small. This applies for symmetric potentials, such as the square and 
parabolic quantum wells. However, for asymmetric quantum wells, the average electric field is non-zero in the direction 
perpendicular to the plane of the 2DEG and is called the interface or quantum well electric field. For experimentally 
achievable semiconductor heterostructures, this field can be as high as 10^ V/cm. Therefore, from Eq. there should 
be an additional (compared with the infinite 3D crystal) mechanism of spin-orbit coupling associated with this field 
and is usually referred to as the Rashba SOI for quantum wells.— When we take into account that the quantum well 
electric field is perpendicular to the heterojunction interface, the spin-orbit Hamiltonian has a contribution which can 
be written for the Rashba coupling as 

n'so = ^ (f^ X P)^ (3) 

within the zero z-component (stationary situation, no electron transfer across the interface). The constant or in 
Eq. which will be simply denoted as a thereafter in this paper, includes universal constants from Eq. ^ and it is 



^so = 



X p 
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proportional to the the interface electric field. The value of a determines the contribution of the Rashba spin-orbit 
coupling to the total electron Hamiltonian. This constant may have values running from (1 — 10) meV- nm. 

Within the single-band effective mass approximation— i^S, the total Hamiltonian of a quasi-one-dimensional electron 
system (QIDES) can be written as 



^ = ^ + ^c(r)+?iso (4) 

where the electron effective mass m* incorporates both the crystal lattice and interaction effects. The form of 
the Hamiltonian derived from the relativistic 4x4 Dirac equation is similar to that which follows from the 8x8 
k • p Hamiltonian^^. Moroz and Barnes^ chose the lateral confining potential T4(r) as a parabola which would be 
appropriate for very narrow wires since the electrons would be concentrated at the bottom of the potential. Such 
narrow QIDES are difficult to achieve experimentally. We are not aware of any experimental evidence or measurement 
of the features arising from the spin-orbit coupling resulting from the parabolic confining potential employed by Moroz 
and Barnes^. So, in this paper, we explore the effects of lateral confinement in which the electrons are essentially free 
over a wide range except close to the edges where the potential rises sharply to confine them. The in-plane electric 
field Ec(r) associated with Vc{v) is given by Ec(r) = —^Vc{v). We assume that the SOI Hamiltonian in Eq. ^ is 

formed by two contributions: T-Lso — ^so' + ^SC)- "^^^ '^'"'^^ "^^^j ^SC) ' ®] a-i'ises from the asymmetry of the 

quantum well, i.e., from the Rashba mechanism for the spin-orbit coupling. For convenience, in what follows we 
will refer to the Rashba mechanism of the spin-orbit coupling as a-coupling. If the lateral confinement is sufficiently 
strong, for narrow and deep potentials or sharp and high potentialsat the edges, then the electric field associated with 
it may not be negligible compared with the interface-induced (Rashba) field. We use 



V.,.) = ,,{„fc(^)+e*(^^)} (5) 

for a conducting channel of width W with well depth Vq. Here, erfc(a;) is the complimentary error function. Plots of 
T4(x)/Vbas a function of x/W are shown in Fig. [2] for three values of W /Iq. For this potential, the Hamiltonian ^ 
gives a term 



where each Gaussian has width £o at the edges a: = and x — W. In Eq. ([6]), J-{x) is related to the electric field 
due to lateral confinement in the x direction. Since ^ W characterizes the steepness of the potentials at the two 
edges, we are at liberty to use a range of values of the ratio of these two lengths, keeping in mind that the in-plane 
confinement must be appreciable if the /3-term is to play a role. Therefore, in most of our calculations, we use only 
one small value of ^o/W to illustrate the effects arising from our model on the conductance and thermoelectric power. 
We introduced the parameter /3o — h^Vo/ {4:\/2tt m* ^c^W), which is expressed in terms of fundamental constants as 
well as Vq and W. The /3o is another Rashba parameter due to the electric confinement along the x direction, and it 
is simply denoted as /3 thereafter in this paper. Comparison of typical electric fields originating from the quantum 
well and lateral confining potentials allows one to conclude that a reasonable estimate for f3 should be roughly 10% 
of a. The /3-SOI term in Eq. ([6]) is asymmetric about the mid-plane a; = W/2 and varies quadratically with the 
displacement from either edge. In this quasi-square well potential, the electron wave functions slightly penetrate into 
the barrier regions. However, we only need energy levels for the calculations of ballistic transport electrons, not the 
wave functions, if we assume electronic system is a spatially-uniform quasi-one-dimensional one. 

The eigenfunctions for the nanowire have the form 
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Since the nanowire is translationally invariant in the y-direction with ky ~ {2t:/ Ly)n, where Ly is a normalization 
length and n = 0, ±1, ±2, • • •, we must solve for ^Pa{x) and tpsix) in Eq. ([7]) numerically due to the presence of edges 
at a; = and x = W . Substituting the wave function in Eq. ([7]) into the Schrodinger equation, i.e. 'H(p{r) = eip{r) 
with e being the eigenenergy, we obtain the two coupled equations 
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2m* \dx'^ ^ ) \dx 
f d^ , - . f d 



k ipA{x) +a\ — + ky\ i})b{x) - likyF{x) tpA{x) = eipAix) , 



2m* \dx'^ ^''^^^^ ~ ^ \dx ~ + l^h^i^)'^B{x) = eipsix) . (8) 

In the absence of any edges, we may simply set J^{x) = for a quantum well, and get ^Pa{x) — Ae^''^^ and 
tpsix) — /Be**^^^, where A and B are independent of x, and k^ is the electron wave number along the x direction, 
which then yields a pair of simultaneous algebraic equations for states A and B. But, in the case when there exist 
edges, we have a pair of coupled differential equations to solve for ipA and ipB which may be analyzed when only /3 is 
not zero and then when both Rashba parameters are non-zero. 

Two parameters of interest are 

£a = h^l2m*a , = h^/2m*l3 , (9) 

with three ratios 

- W/^a , Tl3 = W/£p , To = W/£o . (10) 

In our numerical calculations below, we will use three ratios to determine how sharp the nanowire potential is and 
how strong the Rashba parameters are. 

A. Energy Bands for a — 

When we set a = in Eq. ([5]), iPa{x) and ipsix) are equal to each other {ky —ky) and are solutions of a 
Schrodinger equation with a potential term present. When we solved the wave equations, we imposed the condition 
that the wave functions must vanish when either x or x ^ W holds. However, our calculations showed that the 
wave functions are negligible on the two edges of the nanowire when the confining potential is deep and sharp. The 
effect of the potential depends on ky, leading to a dependence of the transverse energy e^, = e — h^ky/2m* on the 
longitudinal wave number ky. If we set a equal to zero in Eq. ([S]), the solutions are approximately those for a quasi- 
square well when £q/W is chosen small and the potential barriers are high (see Fig. [2]). In this case, the transverse 
energy eigenvalues are approximately given by En = n^Eo, where n — 1, 2, ■ ■ ■ and Eq = n^h^ / {2m*yV^) . In Fig. [31 
we present the energy bands for the calculated transverse energy Ex — e — h'^ky/2m* in units of Eq when the electric 
field-induced Rashba SOI parameter a is set equal to zero so that only the effect from the /3-term is included. The 
two equations coincide and, of course, there is no effect from the Rashba term on Ex when ky = 0, and the eigenvalues 
are equal to those of the square well with high barriers corresponding to io/W ^ 1- However, as kyW is increased. 
Ex decreases linearly as a function of ky before its first drop. On the other hand, the dependence of the transverse 
energy has been shown to be a strong non-linear function of ky at its drop. As a matter of fact, the levels anti-cross 
at a value of kyW which is determined by the chosen value of and io/W. 

The displacing effect of the /3-coupling on the eigenstates is reminiscent of the role played by magnetic field on the 
eigenenergies of a quasi-one-dimensional electron gas with harmonic confinement. However, the anti-crossing seems 
to be a unique property of the square barrier model with high potentials at the edges since it was not reported by 
Moroz and Barnes '* for a parabolic confinement. Figure [3] shows that when (q/W ^ 1, it does not matter what 
value is chosen because the energy level dependence on ky remains the same. When kyW ^ 1, the influence on the 
ky dispersion can be seen. However, this part of the energy spectrum makes a negligible contribution to the transport 
and thermoelectric power. 

B. Energy Bands for a ^ 

In Fig. 131 we present the transverse energies Ex in units of Eq as functions of kyW for a non-zero value of a. The 
plots compare the results for two chosen values of = (black curve) and = 10 (red curve). The energy bands 
for the infinite 2DEG at fc^; = due only to the a-coupling (without transverse confinement) consist of a pair of 
spin-split upward-curved parabolic-like energy dispersions which are displaced in fcj,-space and degenerate at ky = 0, 
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and become split as ky increases. In the presence of edges for a nanowire, there is a discrete set of eigenstates for 
a chosen value of ky. These energy subbands then anti-cross. This anti-crossing effect increases as the value of /3 is 
increased, as demonstrated by comparing the results of Fig.lH These results are qualitatively in agreement with Fig. 3 
in Ref. (see also Ref. [g^) The difference is that for reasonable values of Tp and Tq, the scaled ky becomes larger in 
comparison with a parabolic confinement for one to see the anti-crossing behavior. 

Figures [5] show plots of total energy levels e with tq = 10^ as functions of kyW for Tq, = 5.0, = 1.0 (upper-left 
panel) and — 10.0, Tp — 1.0 (upper-right panel), respectively. For the sake of comparison, the plot for Tq = 0, 
Tp = 2.0, and tq — 10'^ (lower panel) is also presented in this figure. From Figs. [51 we easily find that the energy 
dispersion is symmetric with respect to ky = 0. Our results (two upper panels) further show that, as expected from 
Eq. ([5]) , the a-term has an effect when ky ^ Q 0T^ the unperturbed energy eigenvalues in the absence of any Rashba 
SOI. As ky is increased, the SOI lifts the degeneracy of the energy subbands, as shown by a pair of solid and dashed 
curves degenerated at fcj, = in the upper two panels. We further observed that when the a-coupling is weak (upper- 
left panel), each branch of the energy curves, except for the first branch, has only one local minimum and this occurs 
at ky — Q. As a is increased (upper-right panel), two local minima develop symmetrically on either side oi ky — Q 
with a local maximum at fcj, = for each branch. 



III. MODEL FOR BALLISTIC CHARGE TRANSPORT 



In this section, we briefly outline our method of calculation for the ballistic conductance and electron-diffusion 
thermoelectric power. In our case, the energy bands are symmetric with respect to the wave numbers ±fcy. For the 
symmetric bands, energy dispersion Sj^ky , the Fermi function /o {sj^ky ) , and the group velocity Vj^ky satisfy the relations: 
^j,ky = Sj-ky, h{£j,ky) = foisj-ky), aud Vj^ky = -fc„- Thcrcfore, one can write the following equation in a 
form which includes only positive values of the wave number for the ballistic heat (Q^^-*) and charge (Q^^-*) currents, 
i.e., 



^ eVb(-e)^ ^ 



dfoiej.ky) 



d£j,ky 



de 



j^ky 



(11) 



where ^! = 0, 1, Vb is the bias voltage between the source and drain electrodes, sgn(a;) is the sign function, fco = 0, 
and is the chemical potential. In Eq. (jlip , the whole energy integration performed over the range < fcj, < oo is 
divided into the sum of many sub-integrations between two successive extremum points Sj^k^ for < n < iV, and 
£j,kK is the last minimum point. For each sub-integration over ky, Sj^ky is a monotonic function. In addition, each 
sub-integration in Eq. ([TT]) can be calculated analytically, leading to the following expression for electron-diffusion 
thermoelectric power 



Q(i) fcB 



rQ(0) eg 



(12) 



where T is the temperature, /3 = l/fceT, and the dimensionless conductance g is given by 



.9 = 51 Q,«/ofe,fcJ ■ (13) 

Physically, the quantity g defined in Eg. (|13p represents the number of pairs of the Fermi points at T = OK. In 
Eqs. ([T^ and ([T^ . the summations over n are for all the energy-extremum points on each jth spin-split subband in 
the range < fcj, < oo. The quantity £j,fc„ is the energy at the extremum point ky = For a given jth spin-split 
subband, Cj_„ = 1 (or Cj^„ — —1) for a local energy minimum (maximum) point. The physical conductance G is 
related to g for spin-split subbands through 



(14) 
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IV. NUMERICAL RESULTS OF CHARGE TRANSPORT 

In Figs.m we have displayed comparisons of modified electron density (niD) dependence of the ballistic conductance 
(G) and the electron-diffusion thermoelectric power (S'd) by the a-term in the SOI when T = 4K and W = 568.7 A. 
From the upper panel of Fig.lBl we find that, as a = (black curve), a number of steps in G show up as a result of 
successive populations of more and more spin-degenerate subbands, shown in the lower panel of Fig. [5] In addition, 
the observed plateau becomes wider and wider as higher and higher subbands are occupied by electrons due to 
increased energy-level separation, resulting from the high potential barriers at the two edges. The finite-temperature 
effect can easily be seen from the smoothed steps in this figure. As a is increased to 0.5eV-A (red curve), the steps 
are rightward shifted to higher electron densities due to an enhanced density-of-states from the flattened subband 
dispersion curves by SOI, as seen from the upper-left panel of Fig.O However, the step sharpness remains constant 
in this case. Furthermore, there exists no "pole-like feature" — in this figure, which can be traced back to the absence 
of spike-like feature in the subband dispersion curves, leading to additional local energy minimum/maximum points. 
The suppressed spike-like feature in the subband dispersion curves can be explained by a nonlinear x dependence near 
the center {x = yV/2) of a transverse symmetric potential well with a large /3 value in our model for wide quantum 
wires, instead of a linear x dependence close to the center of the confining potential in the model proposed by Moroz 
and Barnes— for narrow quantum wires. We also see sharp peaks in Sd from the lower panel of Fig. [5] as a = (black 
curve), corresponding to the steps in G, which again comes from successive population of spin-degenerate subbands 
with increased nu^.^ The center of a plateau in G aligns with the minimum of S'd between two peaks. The peaks 
(red curve) are rightward shifted accordingly in electron density when a finite value of a is assumed. 

We have compared in Figs. [7] the results of G and Sd for two values of wire width W at T = 4K and a = 0.5 eV- A. 
We find from the upper panel of Fig.[7]that, as W decreases from 1137.4 A (black curve) to 568.7 A (red curves), the 
steps in G are leftward shifted in electron density, and meanwhile, the steps become sharpened. The step shifting is 
a result of the reduction of SOI effect due to a smaller value for Tq (proportional to aW) with a fixed value of a, 
by comparing the upper-right panel with the upper-left panel of Fig. [5] This leads to a leftward shift in steps for the 
same reason given for the upper panel of Fig. [51 The step sharpening, on the other hand, comes from the significantly 
increased subband separation (proportional to l/W^), which effectively suppresses the thermal-population effect on 
G for smoothing out the conductance steps. The shifting of steps in G with W is also reflected in Sa, as shown in the 
lower panel of Fig.[71 The peaks of S'd get sharpened due to the suppression of Sd in the density region corresponding 
to the widened plateaus of G. 

In order to achieve an overview for the variations of G and Sd with both T and tt-id, we present two contour plots 
of these quantities, respectively, in Figs. [8] with W — 1137. 4A and a — 0.5eV-A. From the upper-left panel of Fig.|8l 
we find that G decreases with T, but increases with riiD in general. The increase of G with riiD is a direct result 
of opening more conduction channels, i.e. more populated subbands, for ballistically-transported electrons. The 
reduction of G with T is a consequence of the dramatic decrease of the chemical potential with T for a fixed tt-id (not 
shown), and then, the decrease of the Fermi function in Eq. (fT3|) for T > 10 K. However, G does increase with T at 
riiD = 4 X lO^cm"^ within the range T < 10 K, as shown in the lower-left panel of Fig. [51 because of the anomalous 
increase of the chemical potential with T in this temperature range whenever the Fermi energy at T = K is set close 
to a minimum of any one of spin-split subbands. When tq is reduced from 10'^ (red curve) to 10 (blue curve) for softer 
potential edges, G is only slightly decreased at low T, but is significantly increased at high T. We also find, from the 
upper- right panel of Fig.[Sl that Sd increases with T but decreases with tt-id. The decrease of Sd with tiyd is simply 
due to the increase of G. As expected, Sd varies linearly with T for small values of T but deviates from the linear 
behavior for large values of T, as shown in the lower-right panel of Fig.|51 However, Sd, in this case, goes towards 
its minimum value about zero at niD = 4 x 10^ cni^^ as T — > because G approaches the quantized value 6 /h for 
this electron density. When tq is reduced from lO'^ (red curve) to 10 (blue curve), Sd is enhanced at high T. 

Finally, we show another pair of contour plots for G and Sd in Fig. [S] as functions of T and niD with W = 568.7 A 
and a = 0.5eV-A. Here, similar to the two upper panels of Fig.[8l G (upper-left panel) increases with riiD but with 
much steeper steps, and Sd (upper-right panel) decreases with niD but at a more rapid rate at a higher temperature. 
At a higher electron density riiD = 10.5 x 10^ cm~^ in the lower-left panel of Fig. [HI for a reduced value of W, the 
range for anomalous increase of G with T expands up to 40 K. This is further accompanied by a locking of G to its 
quantized value at 8 e^/ft. for T < 10 K. However, G is reduced at high T in this case when tq is reduced from 10'^ (red 
curve) to 10 (blue curve). Moreover, the locked G value leads to an almost zero value for Sd within this temperature 
range, as seen from the lower-right panel of Fig. [HI A clear linear dependence of Sd on T is found for T > 40 K due to 
suppressed thermal-population effect by a reduced wire width. When tq is reduced from 10"^ (red curve) to 10 (blue 
curve), no significant changes to Sd can be observed here. 
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V. CONCLUDING REMARKS 



In this paper we investigated the effect that tlie spin-orbit interaction has on the energy band structure, the con- 
ductance and the electron-diffusion thermoelectric power of a nanowire. We used a model in which edge effects for the 
nanowires are taken into account by solving numerically Dirac's equation in a quasi-square potential. Comparing our 
model with the already published work where harmonic confinement was used to describe the transverse confinement, 
we found that the energy bands are different and in addition to crossing effect of the transverse energy bands, there is 
also anticrossing for specific finite values of kyW. The /3-term of the Hamiltonian causes a displacement and a defor- 
mation of the transverse energy band structure which is more pronounced for large values of kyW. The conductance 
plateaus become wider when the electron density is increased as a result of larger energy-level separation as the higher 
subbands are occupied by electrons. Also, due to the nonlinear dependence of /3-SOI on position close to the well 
center, there is no "pole-like feature" in G which is obtained when a harmonic potential is assumed for confinement. 
The electron-diffusion thermoelectric power Sd displays a peak whenever a spin-split subband is populated. At low 
temperature, the variation of Sd is linear in T but deviates from the linear behavior for large values of T. 

We note that if the effect due to electron interaction is strong in our system, the ballistic model in our paper 
cannot be justified. For high-mobility semiconductor quantum wires, electron transport is expected to be ballistic if 
the wire length is shorter than the mean free path of electrons in the system. On the other hand, electron transport 
can also be diffusive if the wire length is longer than the mean free path of electrons but less than the localization 
length. In the latter case, the interaction of electrons with impurities, roughness, phonons and other electrons will 
play a significant role in both the temperature and density dependence of electron mobility in quantum wires. For 
high temperature and low density, electron-phonon scattering is dominant. However, the electron-electron scattering 
becomes significant in a system with high density at low temperature. In this paper, we restricted ourselves to the 
ballistic regime for a short quantum wire. 

Furthermore, In the absence of electron-electron interaction, an effect of lateral spin-orbit coupling on transport is 
to trigger a spontaneous but negligible spin polarization in the nanowire. However, the spin polarization is enhanced 
substantially when the effect of electron-electron interaction is included. The spin polarization may be strong enough 
to result in the appearance of a conductance plateau at a fractional value ('-^ 0.7) of 2e^//i in the absence of any 
external magnetic fieldi^'^ '^'^i^^ The role played by electron-electron interaction on our results for the thermoelectric 
power in a nanowire may also be investigated using finite-temperature Green's function techniques or field theoretic 
methods which would allow the classification of contributing diagrams.— This study should shed some light on whether 
electron-electron interaction enhances the thermoelectric power and how the SOI underlies the Peltier effect, i.e., the 
flow of entropy current in addition to the familiar charge current in an electric field. The effect of electron-electron 
interaction is expected to be significant at very low temperatures (less than IK) and high electron densities. However, 
our model calculations are performed at temperatures at or higher than 4K. 
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FIG. 1: (Color online) Schematic illustration of the nanowire of 2DEG between a source and drain. 
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FIG. 2: (Color online) Plots of Vc{x)/Vo, defined in Eq. (O, as a function of x/W for W/lo = 10 (green curve), W/lo = 50 
(blue curve), and W/lo = 100 (red curve). 
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FIG. 3: (Color online) The transverse energy levels Ex scaled with respect to -Bo, the ground state energy in the absence of any 
SOI, as a function of fcyW for = 0.1 and Tq = 0. 




FIG. 4: (Color online) The transverse energy levels Sx scaled with respect to Eo as a function of fcyW with Tq = 0.3 and 
To = 10"^ for Tp =Q (black curve) and rg = 1.5 (red curve). 
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FIG. 5: (Color online) Plots of total energy levels e as functions of scaled wave number, kyW, for Tc = 5.0 and rp = 1.0 
(upper-left panel) and Tq — 10.0 and — 1.0 (upper-right panel), respectively. Solid and dashed curves represent a pair of 
spin-split subbands in a branch. Here, the plot of spin-degenerated energy levels e for Tc = 0, t/j = 2.0 (lower panel) is also 
included for the comparison. In our calculations, to = 10^ is chosen for these three figures. 
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FIG. 6: (Color online) Comparisons of the conductance G (upper panel), as well as the electron-difltusion thermoelectric power 
Sd (lower panel), as a function of the electron density nio with a wire width W = 568.7 A and a temperature T = 4K for 
Of = (black curves) and a = 0.5eV-A (red curves), respectively. Here, = 1.0 and ro = 10"^ are chosen for the calculations 
in these two figures. 
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FIG. 7: (Color online) Comparisons of G (upper panel) and Sd (lower panel) as a function of niD at T = 4K and with 
a = 0.5eV-A for W = 1137.4 A (black curves) and W = 568.7 A (red curves), separately. Here, = 1.0 and ro = 10^ are 
chosen for the calculations in these two figures. 
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FIG. 8: (Color online) Contour plots of G (upper left) and Sd (upper right) as functions of both T and niD with a = 0.5eV-A 
and W = 1137.4 A. Here, ro — 10"^ is chosen for the upper two Contour plots. The plots of G (lower left) and Sd (lower right) 
are also shown in this figure as a function of T for niD = 4 x 10® cm~^ with tq — 10^ (red curves) and ro = 10 (blue curves). 
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FIG. 9: (Color online) Contour plots of G (upper left) and Sd (upper right) as functions of T and niD with a = 0.5 eV-A and 
W = 568.7 A. Here, tq = 10"^ is chosen for the upper two Contour plots. For a fixed value of niD ~ 10.5 x 10® cm~^, the plots 
for G (lower left) and Sd (lower right) are also shown in this figure as a function of T with to — 10^ (red curves) and tq = 10 
(blue curves). 



